Stereochemical analysis of the antigenic tip of the V3 loop peptide of HIV-1 gp120

A novel multiple turn conformation has been observed for a segment GPGRAFY in the crystal structure of a complex of HIV-1 gp120 V3 loop peptide with the Fab fragment of a neutralizing antibody [Ghiara ct al. (1994) Science 264, 82-85]. A structural motif has been defined for the peptide segment, employing idealized backbone conformations characterized by ranges of virtual C-alpha torsion angles and bond angles. A search of 122 high-resolution protein crystal structures has permitted identification of 24 examples of similar structural motifs. Two major conformational families have been identified, which differ primarily in the conformation at residue 3. The observed conformation at residue 3 in family 1 is left-handed helical (alpha(L)) and that in family 2 is right-handed helical (alpha(R)). Of the 10 examples in family 1, 9 examples have Gly residues at position 3. Of the 12 examples in family 2, 7 examples have Asn/Asp at position 3. Computer modeling of the V3 loop tip sequence using the two backbone conformational families as starting points leads to minimum-energy conformations in which antigenically important side-chains occupy similar spatial arrangements. This stereochemical analysis of the V3 loop tip sequence suggests a rational basis for the design of synthetic analog peptides for use as viral antagonists or synthetic antigens. (C) Munksgaard 1995.

Design of a Non-glycosylated Outer Domain-derived HIV-1 gp120 Immunogen That Binds to CD4 and Induces Neutralizing Antibodies

The outer domain (OD) of the HIV-1 envelope glycoprotein gp120 is an important target for vaccine design as it contains a number of conserved epitopes, including a large fraction of the CD4 binding site.Attempts to design OD-based immunogens in the past have met with little success. We report the design and characterization of an Escherichia coli-expressed OD-based immunogen (ODEC), based on the sequence of the HxBc2 strain. The ODEC-designed immunogen lacks the variable loops V1V2 and V3 and incorporates 11 designed mutations at the interface of the inner and the outer domains of gp120. Biophysical studies showed that ODEC is folded and protease-resistant, whereas ODEC lacking the designed mutations is highly aggregation-prone. In contrast to previously characterized OD constructs, ODEC bound CD4 and the broadly neutralizing antibody b12 but not the non-neutralizing antibodies b6 and F105. Upon immunization in rabbits, ODEC was highly immunogenic,and the sera showed measurable neutralization for four subtype B and one subtype C virus including two b12-resistant viruses. In contrast,sera from rabbits immunized with gp120 did not neutralize any of the viruses. ODEC is the first example of a gp120 fragment-based immunogen that yields significant neutralizing antibodies.

Estimating Frequencies of Minority Nevirapine-Resistant Strains in Chronically HIV-1-Infected Individuals Naive to Nevirapine by Using Stochastic Simulations and a Mathematical Model

Nevirapine forms the mainstay of our efforts to curtail the pediatric AIDS epidemic through prevention of mother-to-child transmission of HIV-1. A key limitation, however, is the rapid selection of HIV-1 strains resistant to nevirapine following the administration of a single dose. This rapid selection of resistance suggests that nevirapine-resistant strains preexist in HIV-1 patients and may adversely affect outcomes of treatment. The frequencies of nevirapine-resistant strains in vivo, however, remain poorly estimated, possibly because they exist as a minority below current assay detection limits. Here, we employ stochastic simulations and a mathematical model to estimate the frequencies of strains carrying different combinations of the common nevirapine resistance mutations K103N, V106A, Y181C, Y188C, and G190A in chronically infected HIV-1 patients naive to nevirapine. We estimate the relative fitness of mutant strains from an independent analysis of previous competitive growth assays. We predict that single mutants are likely to preexist in patients at frequencies (similar to 0.01% to 0.001%) near or below current assay detection limits (>0.01%), emphasizing the need for more-sensitive assays. The existence of double mutants is subject to large stochastic variations. Triple and higher mutants are predicted not to exist. Our estimates are robust to variations in the recombination rate, cellular superinfection frequency, and the effective population size. Thus, with 10(7) to 10(8) infected cells in HIV-1 patients, even when undetected, nevirapine-resistant genomes may exist in substantial numbers and compromise efforts to prevent mother-to-child transmission of HIV-1, accelerate the failure of subsequent antiretroviral treatments, and facilitate the transmission of drug resistance.

Taking Multiple Infections of Cells and Recombination into Account Leads to Small Within-Host Effective-Population-Size Estimates of HIV-1

Whether HIV-1 evolution in infected individuals is dominated by deterministic or stochastic effects remains unclear because current estimates of the effective population size of HIV-1 in vivo, N-e, are widely varying. Models assuming HIV-1 evolution to be neutral estimate N-e similar to 10(2)-10(4), smaller than the inverse mutation rate of HIV-1 (similar to 10(5)), implying the predominance of stochastic forces. In contrast, a model that includes selection estimates N-e>10(5), suggesting that deterministic forces would hold sway. The consequent uncertainty in the nature of HIV-1 evolution compromises our ability to describe disease progression and outcomes of therapy. We perform detailed bit-string simulations of viral evolution that consider large genome lengths and incorporate the key evolutionary processes underlying the genomic diversification of HIV-1 in infected individuals, namely, mutation, multiple infections of cells, recombination, selection, and epistatic interactions between multiple loci. Our simulations describe quantitatively the evolution of HIV-1 diversity and divergence in patients. From comparisons of our simulations with patient data, we estimate N-e similar to 10(3)-10(4), implying predominantly stochastic evolution. Interestingly, we find that N-e and the viral generation time are correlated with the disease progression time, presenting a route to a priori prediction of disease progression in patients. Further, we show that the previous estimate of N-e>10(5) reduces as the frequencies of multiple infections of cells and recombination assumed increase. Our simulations with N-e similar to 10(3)-10(4) may be employed to estimate markers of disease progression and outcomes of therapy that depend on the evolution of viral diversity and divergence.

Estimating the Threshold Surface Density of Gp120-CCR5 Complexes Necessary for HIV-1 Envelope-Mediated Cell-Cell Fusion

Reduced expression of CCR5 on target CD4(+) cells lowers their susceptibility to infection by R5-tropic HIV-1, potentially preventing transmission of infection and delaying disease progression. Binding of the HIV-1 envelope (Env) protein gp120 with CCR5 is essential for the entry of R5 viruses into target cells. The threshold surface density of gp120-CCR5 complexes that enables HIV-1 entry remains poorly estimated. We constructed a mathematical model that mimics Env-mediated cell-cell fusion assays, where target CD4(+)CCR5(+) cells are exposed to effector cells expressing Env in the presence of a coreceptor antagonist and the fraction of target cells fused with effector cells is measured. Our model employs a reaction network-based approach to describe protein interactions that precede viral entry coupled with the ternary complex model to quantify the allosteric interactions of the coreceptor antagonist and predicts the fraction of target cells fused. By fitting model predictions to published data of cell-cell fusion in the presence of the CCR5 antagonist vicriviroc, we estimated the threshold surface density of gp120-CCR5 complexes for cell-cell fusion as similar to 20 mu m(-2). Model predictions with this threshold captured data from independent cell-cell fusion assays in the presence of vicriviroc and rapamycin, a drug that modulates CCR5 expression, as well as assays in the presence of maraviroc, another CCR5 antagonist, using sixteen different Env clones derived from transmitted or early founder viruses. Our estimate of the threshold surface density of gp120-CCR5 complexes necessary for HIV-1 entry thus appears robust and may have implications for optimizing treatment with coreceptor antagonists, understanding the non-pathogenic infection of non-human primates, and designing vaccines that suppress the availability of target CD4(+)CCR5(+) cells.

Designed Cyclic Permutants of HIV-1 gp120: Implications for Envelope Trimer Structure and Immunogen Design

Most HIV-1 broadly neutralizing antibodies are directed against the gp120 subunit of the env surface protein. Native env consists of a trimer of gp120-gp41 heterodimers, and in contrast to monomeric gp120, preferentially binds CD4 binding site (CD4bs)-directed neutralizing antibodies over non-neutralizing ones. Some cryo-electron tomography studies have suggested that the V1V2 loop regions of gp120 are located close to the trimer interface. We have therefore designed cyclically permuted variants of gp120 with and without the h-CMP and SUMO2a trimerization domains inserted into the V1V2 loop. h-CMP-V1cyc is one such variant in which residues 153 and 142 are the N- and C-terminal residues, respectively, of cyclically permuted gp120 and h-CMP is fused to the N-terminus. This molecule forms a trimer under native conditions and binds CD4 and the neutralizing CD4bs antibodies b12 with significantly higher affinity than wild-type gp120. It binds non-neutralizing CD4bs antibody F105 with lower affinity than gp120. A similar derivative, h-CMP-V1cycl, bound the V1V2 loop-directed broadly neutralizing antibodies PG9 and PG16 with similar to 20-fold higher affinity than wild-type JRCSF gp120. These cyclic permutants of gp120 are properly folded and are potential immunogens. The data also support env models in which the V1V2 loops are proximal to the trimer interface.

Molecular Epidemiology of HIV-1 Subtypes in India: Origin and Evolutionary History of the Predominant Subtype C

Background: India has the third largest HIV-1 epidemic with 2.4 million infected individuals. Molecular epidemiological analysis has identified the predominance of HIV-1 subtype C (HIV-1C). However, the previous reports have been limited by sample size, and uneven geographical distribution. The introduction of HIV-1C in India remains uncertain due to this lack of structured studies. To fill the gap, we characterised the distribution pattern of HIV-1 subtypes in India based on data collection from nationwide clinical cohorts between 2007 and 2011. We also reconstructed the time to the most recent common ancestor (tMRCA) of the predominant HIV-1C strains. Methodology/Principal Findings: Blood samples were collected from 168 HIV-1 seropositive subjects from 7 different states. HIV-1 subtypes were determined using two or three genes, gag, pol, and env using several methods. Bayesian coalescent-based approach was used to reconstruct the time of introduction and population growth patterns of the Indian HIV-1C. For the first time, a high prevalence (10%) of unique recombinant forms (BC and A1C) was observed when two or three genes were used instead of one gene (p<0.01; p = 0.02, respectively). The tMRCA of Indian HIV-1C was estimated using the three viral genes, ranged from 1967 (gag) to 1974 (env). Pol-gene analysis was considered to provide the most reliable estimate 1971, (95% CI: 1965-1976)]. The population growth pattern revealed an initial slow growth phase in the mid-1970s, an exponential phase through the 1980s, and a stationary phase since the early 1990s. Conclusions/Significance: The Indian HIV-1C epidemic originated around 40 years ago from a single or few genetically related African lineages, and since then largely evolved independently. The effective population size in the country has been broadly stable since the 1990s. The evolving viral epidemic, as indicated by the increase of recombinant strains, warrants a need for continued molecular surveillance to guide efficient disease intervention strategies.

Stochastic Simulations Suggest that HIV-1 Survives Close to Its Error Threshold

The use of mutagenic drugs to drive HIV-1 past its error threshold presents a novel intervention strategy, as suggested by the quasispecies theory, that may be less susceptible to failure via viral mutation-induced emergence of drug resistance than current strategies. The error threshold of HIV-1, mu(c), however, is not known. Application of the quasispecies theory to determine mu(c) poses significant challenges: Whereas the quasispecies theory considers the asexual reproduction of an infinitely large population of haploid individuals, HIV-1 is diploid, undergoes recombination, and is estimated to have a small effective population size in vivo. We performed population genetics-based stochastic simulations of the within-host evolution of HIV-1 and estimated the structure of the HIV-1 quasispecies and mu(c). We found that with small mutation rates, the quasispecies was dominated by genomes with few mutations. Upon increasing the mutation rate, a sharp error catastrophe occurred where the quasispecies became delocalized in sequence space. Using parameter values that quantitatively captured data of viral diversification in HIV-1 patients, we estimated mu(c) to be 7 x 10(-5) -1 x 10(-4) substitutions/site/replication, similar to 2-6 fold higher than the natural mutation rate of HIV-1, suggesting that HIV-1 survives close to its error threshold and may be readily susceptible to mutagenic drugs. The latter estimate was weakly dependent on the within-host effective population size of HIV-1. With large population sizes and in the absence of recombination, our simulations converged to the quasispecies theory, bridging the gap between quasispecies theory and population genetics-based approaches to describing HIV-1 evolution. Further, mu(c) increased with the recombination rate, rendering HIV-1 less susceptible to error catastrophe, thus elucidating an added benefit of recombination to HIV-1. Our estimate of mu(c) may serve as a quantitative guideline for the use of mutagenic drugs against HIV-1.

Conformational dynamics of HIV-1 protease: a comparative molecular dynamics simulation study with multiple amber force fields

Flap dynamics of HIV-1 protease (HIV-pr) controls the entry of inhibitors and substrates to the active site. Dynamical models from previous simulations are not all consistent with each other and not all are supported by the NMR results. In the present work, the er effect of force field on the dynamics of HIV-pr is investigated by MD simulations using three AMBER force fields ff99, ff99SB, and ff03. The generalized order parameters for amide backbone are calculated from the three force fields and compared with the NMR S2 values. We found that the ff99SB and ff03 force field calculated order parameters agree reasonably well with the NMR S2 values, whereas ff99 calculated values deviate most from the NMR order parameters. Stereochemical geometry of protein models from each force field also agrees well with the remarks from NMR S2 values. However, between ff99SB and ff03, there are several differences, most notably in the loop regions. It is found that these loops are, in general, more flexible in the ff03 force field. This results in a larger active site cavity in the simulation with the ff03 force field. The effect of this difference in computer-aided drug design against flexible receptors is discussed.

Design of an Escherichia coli expressed HIV-1 gp120 fragment Immunogen that binds to b12 and induces broad and potent neutralizing antibodies

b12, one of the few broadly neutralizing antibodies against HIV-1, binds to the CD4 binding site (CD4bs) on the gp120 subunit of HIV-1 Env. Two small fragments of HIV-1 gp120, b121a and b122a, which display about 70% of the b12 epitope and include solubility-enhancing mutations, were designed. Bacterially expressed b121a/b122a were partially folded and could bind b12 but not the CD4bs-directed non-neutralizing antibody b6. Sera from rabbits primed with b121a or b122a protein fragments and boosted with full-length gp120 showed broad neutralizing activity in a TZM-bl assay against a 16-virus panel that included nine Tier 2 and 3 viruses as well as in a five-virus panel previously designed to screen for broad neutralization. Using a mean IC50 cut-off of 50, sera from control rabbits immunized with gp120 alone neutralized only one virus of the 14 non-Tier 1 viruses tested (7%), whereas sera from b121a- and b122a-immunized rabbits neutralized seven (50%) and twelve (86%) viruses, respectively. Serum depletion studies confirmed that neutralization was gp120-directed and that sera from animals immunized with gp120 contained lower amounts of CD4bs-directed antibodies than corresponding sera from animals immunized with b121a/b122a. Competition binding assays with b12 also showed that b121a/2a sera contained significantly higher amounts of antibodies directed toward the CD4 binding site than the gp120 sera. The data demonstrate that it is possible to elicit broadly neutralizing sera against HIV-1 in small animals.

Simulations reveal that the HIV-1 gp120-CD4 complex dissociates via complex pathways and is a potential target of the polyamidoamine (PAMAM) dendrimer

The polyamidoamine (PAMAM) dendrimer prevents HIV-1 entry into target cells in vitro. Its mechanism of action, however, remains unclear and precludes the design of potent dendrimers targeting HIV-1 entry. We employed steered molecular dynamics simulations to examine whether the HIV-1 gp120-CD4 complex is a target of PAMAM. Our simulations mimicked single molecule force spectroscopy studies of the unbinding of the gp120-CD4 complex under the influence of a controlled external force. We found that the complex dissociates via complex pathways and defies the standard classification of adhesion molecules as catch and slip bonds. When the force loading rate was large, the complex behaved as a slip bond, weakening gradually. When the loading rate was small, the complex initially strengthened, akin to a catch bond, but eventually dissociated over shorter separations than with large loading rates. PAMAM docked to gp120 and destabilized the gp120-CD4 complex. The rupture force of the complex was lowered by PAMAM. PAMAM disrupted salt bridges and hydrogen bonds across the gp120-CD4 interface and altered the hydration pattern of the hydrophobic cavity in the interface. In addition, intriguingly, PAMAM suppressed the distinction in the dissociation pathways of the complex between the small and large loading rate regimes. Taken together, our simulations reveal that PAMAM targets the gp120-CD4 complex at two levels: it weakens the complex and also alters its dissociation pathway, potentially inhibiting HIV-1 entry.

Estimating the fraction of progeny virions that must incorporate APOBEC3G for suppression of productive HIV-1 infection

The contest between the host factor APOBEC3G (A3G) and the HIV-1 protein Vif presents an attractive target of intervention. The extent to which the A3G-Vif interaction must be suppressed to tilt the balance in favor of A3G remains unknown. We employed stochastic simulations and mathematical modeling of the within-host dynamics and evolution of HIV-1 to estimate the fraction of progeny virions that must incorporate A3G to render productive infection unsustainable. Using three different approaches, we found consistently that a transition from sustained infection to suppression of productive infection occurred when the latter fraction exceeded similar to 0.8. The transition was triggered by A3G-induced hypermutations that led to premature stop codons compromising viral production and was consistent with driving the basic reproductive number, R-o, below unity. The fraction identified may serve as a quantitative guideline for strategies targeting the A3G-Vif axis. (C) 2013 Elsevier Inc. All rights reserved.

Stabilizing the Native Trimer of HIV-1 Env by Destabilizing the Heterodimeric Interface of the gp41 Postfusion Six-Helix Bundle

The HIV-1 envelope glycoprotein (Env) is a trimer of gp120-gp41 heterodimers and is essential for viral entry. The gp41 subunit in native, prefusion trimeric Env exists in a metastable conformation and attains a stable six-helix bundle (6-HB) conformation comprised of a trimer of N-heptad repeat (NHR) and C-heptad repeat (CHR) heterodimers, that drives the fusion of viral and cellular membranes. We attempted to stabilize native Env trimers by incorporation of mutations at the NHR-CHR interface that disrupt the postfusion 6-HB of gp41. The mutations V570D and I573D stabilize native Env of the HIV-1 JRFL strain and occlude nonneutralizing epitopes to a greater extent than the previously identified I559P mutation that is at the interface of the NHR trimers in the 6-HB. The mutations prevent soluble-CD4 (sCD4)-induced gp120 shedding and 6-HB formation. In the context of cell surface-expressed JRFL Env, introduction of a previously reported additional disulfide between residues A501 and T605 perturbs the native conformation, though this effect is partially alleviated by furin coexpression. The data suggest that positions 570 and 573 are surface proximal in native Env and that the NHR homotrimeric coiled coil in native Env terminates before or close to residue 573. Aspartic acid substitutions at these positions stabilize native trimers through destabilization of the postfusion 6-HB conformation. These mutations can be used to stabilize Env in a DNA vaccine format. IMPORTANCE The major protein on the surface of HIV-1 is the envelope (Env) glycoprotein. Env is a trimer of gp120-gp41 heterodimers. gp120 is involved in receptor/coreceptor binding and gp41 in the fusion of viral and cellular membranes. Like many other viral fusion proteins, the gp41 subunit in native trimeric Env exists in a metastable conformation. gp41 readily forms a stable six-helix bundle (6-HB) conformation comprised of a trimer of N-heptad repeat (NHR) and C-heptad repeat (CHR) heterodimers that drives fusion of viral and cellular membranes. While it is expected that native Env is a good immunogen, its metastability results in exposure of immunodominant nonneutralizing epitopes. In the present study, we stabilize native Env trimers by incorporation of a number of different mutations at the NHR-CHR interface that disrupt the postfusion 6-HB of gp41. The stabilized constructs described here can be incorporated into DNA vaccine candidates.